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Abstract 



We study one and two hole properties of the t-J model at half-fihing on the 
square lattice using series expansion methods at T = 0. The dispersion curve 
for one hole excitations is calculated and found to be qualitatively similar to 
that obtained by other methods, but the bandwidth for small t/J is some 
20% larger than given previously. We also obtain the binding energy and 
dispersion relation for two hole bound states. The lowest bound state as t/J 
increases is found to be first d-wave, and then p-wave, in accordance with 
predictions based upon the Kohn-Luttinger effect. We also make a similar 
study for the t-Jz model. 

PACS Indices: 71.27.+a, 71.10Fd 
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I. INTRODUCTION 



Models of strongly correlated electrons in two dimensions have been much studied in 
recent years as part of a search for a minimal model of the high-Tc superconducting cuprates, 
and other fascinating new antiferromagnetic systems The t-J model, defined by the 

Hamiltonian 

H = -tY. ^(4.9,. + cl,c,,,)P + J^(S, • S, - \n,n,) (1) 

(«i>:0- {ij) 

has become a generic model for mobile electrons/holes in doped antiferromagnets. In this 
expression the c|^, q^o- are the usual electron creation and destruction operators, Sj are 
the electron spin operators and the sum is over nearest neighbors pairs. P is a projection 
operator which ensures that doubly occupied states are excluded. 

The t-J model originally was derived ^ as the strong correlation limit of the Hubbard 
model, with J = 2t'^/U, but it is now generally treated as an effective Hamiltonian in 
its own right, with t and J independent parameters. For the cuprates it is estimated 
that t/J ~ 3, which is well outside the region of validity of the Hubbard mapping. An 
independent "derivation" of the t-J model as an effective Hamiltonian for the cuprates has 
been presented by Zhang and Rice 0. 

At exactly half-filling (one electron per site) the kinetic term vanishes because of the 
single-occupancy restriction and the model reduces to an antiferromagnetic insulator. It is 
known from a variety of numerical and analytic calculations ^ that this model, on the square 
lattice, has nonzero long range antiferromagnetic (Neel) order, reduced by strong quantum 
fluctuations. Removal of a small number of electrons (by doping) will allow mobility of holes 
and will reduce the antiferromagnetic order. The physics of a small number of holes in a 
dynamic antiferromagnetic background remains a challenging problem. Previous studies of 
the t-J model at T = have used a variety of approaches. Exact diagonalization methods 
have been used extensively on clusters up to 32 sites [p|-[TT[], but suffer from substantial finite- 
size corrections and an inability to treat extended excitations. Green's function Monte Carlo 
calculations have been performed for larger lattices, up to 8 x 8 [Q, and show, for example. 



a significant decrease of the two-hole binding energy with increasing lattice size. A variety of 
variational approaches and analytic many body methods have also been employed |[T3|-[T9 . 



The picture that emerges from this body of work is roughly as follows. For one hole the 
ground state energy is given by 

eiJt~-3.2 + 2.8(J/t)°-^3 (2) 

where the coefficients are approximate numerical estimates. The dispersion curve for 1- 
hole states is qualitatively similar from different methods, with minima at (7r/2,7r/2) and 
equivalent points. The bandwidth from the largest diagonalizations scales as 

W/t ~ -0.33 + 2.04( J/t)°-6^ (3) 

where the coefficients are again estimates and are not known with great precision. The 



spectral function for 1-hole states has been calculated |10| and shows a clear quasiparticle 



peak, together with a continuum at higher energies. For two holes the binding energy, 
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spectral function and pair susceptibility have been computed. The results support the 
existence of a bound state with S = and d-wave symmetry for J > J^, with ~ 0.3t 
although not all methods agree and this is still controversial. The dispersion curve for two- 
hole bound states is also of interest — to our knowledge the only calculation of this is due 



to Eder [g. 

The aim of the present work is to study the physics of one-hole and two-hole states 
in the half-filled t-J model at zero temperature via linked-cluster series expansions. This 
approach has been developed extensively in recent years by our group pO|- p^] and others 
Most of the applications to date have been to spin systems, although we have 



included fermions in lattice gauge calculations. Shi and Singh |^ have used this approach 
to study the Hubbard model at half filling, via an expansion about an additional Ising term 
included in the Hamiltonian. For completeness we also mention the use of high-temperature 
expansions to study the t-J model [^,|^. The present work of course addresses different 
questions which cannot be probed by high-temperature series. 

The plan of the paper is as follows. In Sec. II, the series expansion method is briefly 
reviewed, and the results for one-hole properties of the t-J model are presented. In Sec. Ill, 
we study the dispersion relation and the binding energy for two-hole bound states. In Sec. 
IV, we study the one-hole and two-hole properties of the t-J^ model. Sec. V is devoted to 
a summary and discussion. 



II. ONE-HOLE PROPERTIES 

The series expansion method is based on a linked cluster formalism associated with 
standard Rayleigh-Schrodinger perturbation theory. 

The ground state energy Eq for a lattice of sites can be written as a sum over all 
connected clusters {a} as 

E^ = T.C^^a (4) 

a 

where is a geometrical embedding factor and ea is the "cumulant energy" of cluster a. 
The same formula can be used for the energy of any particular cluster /3 

E'o=T.C^e^ (5) 

a 

where the sum is now over all sub-clusters embeddable in (3. The cluster energy is calcu- 
lated perturbatively, then the cumulant energies are obtained iteratively from Eq. and 
substituted into Eq. @ to obtain the final series. The formalism for excited states, needed 
to compute the dispersion curves, is a little more complex P3 . 
We introduce an Ising spin anisotropy 

J(S, ■ S, - ^n,n,) ^ J,iS!S] - ^n,n,) + J,,(5f 5J + SfS^) 

= J^iiS^S; - ^n,n,) + x{S!S^ + SfS])] (6) 

where x = Jxy/Jzi and rewrite the Hamiltonian for the t-J model as 
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HI J, = Ho + xV 



(7) 



with 



^0 = T.iS!S] - -n,n,)+rY.i-iyS! (8) 
+ j:iS!S] + SfSy)-rY.i-iyS! (9) 



with a staggered 2;-field of strength r to improve convergence . This term vanishes in the 
total Hamiltonian in the isotropic exchange hmit x — 1. The Ising part of the exchange, 
Ho, forms the unperturbed Hamiltonian while the perturbation V includes both the hole 
hopping terms and transverse spin fluctuations. The original Hamiltonian for the t-J model 
is recovered in the limit x = 1 (i.e. J = = Jxy)- We expect the series expansion to give 
more reliable results in the small t/ J region, where the perturbation V is less important. 

We directly compute the series for the ground state energy difference between the system 
with one (or two) hole(s) and the system with no hole. Two types of series can be obtained 
in this way. By fixing y = t/J^y and r we can compute series in the single variable x up to 
order L (where L is typically of order 11-13), 

E/J, = j2c,{y,r)x' (10) 

We refer to these as x-series. Extrapolation to x = 1 is carried out by standard Fade 
approximant or integrated differential approximant methods . In the other approach we 
keep both x and t/ Jz as expansion parameters and obtain double series of the form 

E/Jz = Y.d,,{r){t/JzYx^ (11) 

where the coefficients dij are computed up to order i + j < L. Extrapolation to x = 1 then 
yields a series in t/J, 

E/J = Y.e.{t/jy (12) 

i 

which we refer to as a t-series. 

We first consider states with one hole, which may exist either on the up-sublattice (A) or 
the down-sublattice (B), where the up and down sublattice are represented in Fig. |l|(a) by 
open and full circles, respectively. Here (e^;, Cy) are the usual primitive translation vectors. 
The unperturbed wave function for the zero-momentum spin-up state is taken to be 

m = T.^^^\N) (13) 

where |A^) is the classical Neel state. Now a hole in the up (down) sublattice can never 
transform into a single hole in the down (up) sublattice, since the Hamiltonian conserves 
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the total spin J2i Si', so that we should choose {cx', e^/) shown in Fig. [l|(a) as the primitive 
translation vectors, rather than {ex,ey). The unit vectors {ex',ey/) are those of the antifer- 
romagnetic bipartite lattice. The Brillouin zones in reciprocal space corresponding to these 
different choices are shown in Fig. |l|(b). Most previous calculations, however, have used 
{ex,ey) as unit vectors (corresponding to {kx,ky) in momentum space), instead of {ex',ey>), 
so we will also conform to this convention (and set a to one) unless it is specified otherwise. 

For momentum k = we obtain the energy to order L = 13, involving 8739 distinct 
clusters containing up to 13 sites. For k 7^ it is necessary to distinguish between clusters 
having different spatial orientation. We expand the energy to order L = 11, involving 
23546 clusters with up to 11 sites. The series coefficients for y = 0.5 at fixed momentum 
{kx,ky) = (7r/2,7r/2) and different r are listed in Table |l[ where we can see the series 
converges best around r = 2. The extrapolations of these series by using the integrated 
differential approximant methods are shown in Figure]^, where we can see that at a: = 1, the 
extrapolations from different r series agree with each other very well, while the case r ~ 2 
gives the smallest error bar. In Table ||, we list the full spectrum series for y = 0.5 and 
r = 0. The other coefficients are too lengthy to publish but can be supplied on request. 

In Figure ^ we show the energy of a one-hole state for k along symmetry directions in 
the Brillouin zone, for various ratios t/J obtained from x-series. The lowest energy occurs 
at {kx,ky) = (7r/2,7r/2) and equivalent points. This result, and the overall shape of the 



one-hole dispersion curve are in very good agreement with previous work ||13| , |30| , |31| , pl)| , p 
We note that the symmetry and the redundant information evident in Fig. ^ are the results 
of choosing a Brillouin zone which is two times too big. If we choose {kx', kyi) as defining the 
bipartite Brillouin zone, the corresponding dispersion relation along symmetry directions 
in the Brillouin zone is shown in Fig. ^. Figure ^ shows the one-hole bandwidth, W = 
£■(0,0) — i?(7r/2, 7r/2), as a function of J/t, compared with results from other approaches. 
For J/t around 0.5 our results agree well with other work, but for large J/t our result is 
significantly larger than that obtained by Martinez and Horsch ||^. Also we found a peak 
value of W/t ~ 1.15 at J/t ~ 1.25, rather than W/t ~ 1 at J/t = 0.8 as obtained by 
Martinez and Horsch |^, and Liu and Manousakis [|l^ . 

Using the two- variable expansion we obtain the one- hole energy £^(k) at some particular 
values of k: 

E^^^iO, 0)/J = 2.164(2) - 0.655(8)(t/J)^ - 0.031(5)(t/J)^ + 0[{t/J)^] 
E^^(7r/2, 7r/2)/J = 2.164(2) - 3.91(4)(t/J)' + 9.1(4)(t/J)^ + 0[{t/jf] (14) 
E^'^iTT, 0)/J = 2.164(2) - 2.68(2)(t/J)2 + 3.06(4)(t/ J)^ + 0[(t/J)^] 

These coefficients are not exact as they involve an extrapolation to x = 1. The leading 



coefficient 2.164(3) can be compared with 2.1552 from a spin wave calculation and 
2.193(7) from a Monte Carlo simulation P3[ . 



The bandwidth can be obtained as 

W/J = 3.26(5)(t/J)2 - 9.1(4)(t/J)^ + 0[{t/jf] (15) 

This result is also plotted in Figure ^ as a dashed line, and it agrees well at small t/J 
with that obtained from x-series. The leading coefficient 3.26(5) should be compared with 
previous calculations giving a value of 2 |3ll] or 2.6 p4 |. 
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III. TWO-HOLE PROPERTIES 



In the previous section we developed a series expansion approach to calculate one-hole 
properties of the t-J model, and showed how reliable quantitative results could be obtained. 
In the present section we extend this approach to consider states with two holes. 

We start from an unperturbed wave-function consisting of two nearest-neighbour holes 
with opposite spin, and overall momentum k: 

|vl/(k)) = Q,Te^^-''^(Q_e.ae-^'-'^/' + Q+e^ae^'^'^^' + c.-e.^e"^'"'^/' + Q+e„ie''""/') |N) (16) 

Here again there are two different choices of Brillouin zone possible. One may choose to 
use a diagonal unit vector of length aj\pi (where a is the lattice spacing), which connects 
the centre points of two adjacent bonds of the original lattice as shown in Fig. P(a); this 
corresponds to the large, outer Brillouin zone in Fig. ^(b). This was basically the convention 
used by Eder |T^. Alternatively, and more correctly, one may choose a diagonal unit vector 
of length A/2a, differentiating between the A and B sublattices. This gives rise again to 
the Brillouin zone of the bipartite lattice, marked as the inner zone in Fig. ^(b). In this 
second case, there are 4 independent bond states in the unit cell, and one may form linear 
combinations of them making up the conventional s-, pa;-, Vy ci^nd d-wave states. With 
the first choice, all these four states are included within one Brillouin zone: the points 
corresponding to zero-momentum s-, p^^-, Py- and d-wave states are shownQ in Fig. |^(b). 
Then one finds that the different states are related by equations such as 

E'\k^, ky) = E^^'K -k^,7C- ky) (17) 

(in units where a/\/2 has been rescaled to one). This relationship is evident in Fig. 19 of 



Ref. \U 



We have computed series for the two- hole spectrum ii^^'^(k) up to order L = 11. This 
involves 46,440 distinct clusters of up to 12 sites. The x-series coefficients for y = 0.5 and 
r = are listed in Table |T[ 

Using the two- variable expansion, after extrapolation to x = 1, we obtain the following 
estimates for the two- hole s-, p- and c?-wave pair energy: 

E^'^/J = E^\0, 0)/J = 3.798(2) - 7.66(3)(t/ J)^ + 17.5(5)(t/J)^ + 0[(t/J)''] 

E^P/J = E^\0, n)/J = 3.798(2) - 3.705(10)(t/ J)' - 4.1(l)(t/J)^ + 0[{t/jf] (18) 

E^yj = E^\Tt, 7r)/J = 3.798(2) - 0.52(3)(t/J)' - 1.62(5)(t/J)^ + 0[{t/jf] 

So the binding energies defined by 

Eb(k) = E^'^ik) - 2E^'^(7r/2, 7r/2) (19) 

for S-, p- and d-wave pair states are: 



The "pa;-" and "Pj/-" wave states are degenerate by reflection symmetry. The points marked 
and py in Fig. ^(b) are really linear combinations of the proper px and Py states. 
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El'' /J = -0.529(4) + 0.16(10)(t/J)' - 0.7(13)(t/ J)^ + 0{{t/Jf] 

Ely J = -0.529(4) + 4.1(l)(t/J)' - 22.3(9)(t/J)^ + 0[{t/jf] (20) 

Ely J = -0.529(4) + 7.3(l)(t/J)' - 19.8(9)(t/ J)^ + 0[{t/jf] 



The binding energy for t = has previously been given as —0.51 by spin wave theory [32 



and —0.58(2) by Monte Carlo |Q. Our result is in relatively good agreement but suggests 



that the Monte Carlo calculations may overestimate the binding energy somewhat. 

We have also computed the dispersion curves for the two-hole bound state from x-series, 
and these are shown in Figure 0, for various small t/J. The minimum occurs at (7r,7r) 
(d-wave pair), and the maximum occurs at (0, 0) (s-wave pair), while the energy at (0, 0) is 
almost unchanged as t/ J increases. The curves are quite flat for small t/J, as expected, and 
develop more structure as t/J increases. The results are broadly similar to those of Eder 

The S-, p- and rf-wave pair binding energies, obtained from both the x-series and the 
two- variable expansion as a function of t/J, are shown in Figure |^. Note that for x-series, 
the extrapolation is done on series for i?b(k), rather than the 1-hole and 2- hole energies 
individually, because values of these two energies are quite large, while their difference is 
usually small. The convergence is excellent for small t/J, and we believe our results are very 
accurate there. The uncertainty increases for larger t. We can see from this figure that as 
t/J increases from zero, the s-wave binding energy decreases very quickly, and the binding 
disappears around t/J 0.295. The p-wave binding energy also decreases quickly in the 
small t/J region, but it turns out to be very fiat beyond t/J ~ 1, and the binding energy 
is very small. The d-wave binding energy decreases slightly as t/ J increases from zero. Our 
results are consistent with the ci-wave binding disappearing around t/J ~ 2.5. Beyond that 
point, the p-wave state appears to have the lower energy, and remains very weakly bound 
according to our results. 

Figure ^ also shows some other numerical results for the 2-hole binding energy. The solid 
circles with error bars are GFMC results for the 8x8 lattice obtained by Boninsegni and 
Manousakis [|l^], and lie somewhat below our result, though broadly similar in trend: the 



discrepancy may easily be attributed to a finite-size effect. The solid triangles are the results 



of a recent 32-site finite- lattice calculation |3^. The open circles are the results of a linear 



extrapolation in 1/A^ from the finite-lattice data by Chernyshev, Leung and Gooding 
but they discuss reasons why such a naive extrapolation is not to be relied on especially at 
large t/J. The results lie somewhat above our series results for the rf-wave binding energy. 



but again have a broadly similar trend. Chernyshev, Leung and Gooding |^ find that the 
lowest bound state on the 32-site lattice is a singlet state with dx2-y2 symmetry. 

The apparent cross-over between the p-wave and rf-wave bound states came initially as 
a surprise. But in fact it accords very well with predictions [0,0 based on the Kohn- 



Luttinger effect. In a repulsive Fermi liquid, pairing in high-angular momentum states can 
occur due to a singularity in the 2-particle scattering amplitude. This effect was previously 
believed not to occur in two dimensions, because it does not appear up to second order 



in perturbation theory; but it was recently shown [36,37| that a singularity does arise at 



higher order. Chubukov []5^ showed that in the dilute limit the effect survives down to / = 1 
and gives rise to p-wave pairing. He also pointed out that the same effect will occur in an 



arbitrarily dense Fermi liquid. Both numerical [|5H] and analytical pD|,p]| approaches have 
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shown that at low densities j9-wave pairing is dominant for J/t <^ 1, d-wave is dominant 
at J/t > 1, and there is a region of s-wave dominance for J/t > 2. Figure ^ appears to 
indicate that a similar pattern extends even to the region of half-filling. In fact, the earlier 
variational calculation]^ of Eder also shows a crossover between p-wave and d-waye states 
(Fig. 21 of Ref. 0) at around t/J =2. 

Several theoretical approaches [p!7| -[T9|| have predicted that magnon exchange should give 
rise to shallow, long-range bound states even at large t/J values. The behaviour of the 
p-wave state looks consistent with this idea, but not the d-wave state. It is possible that 
another state (or states) in the (i-wave sector could cross over the one we have calculated, 
and become the lowest energy state at large t/J. Such a crossover could typically imply 
a square-root singularity in the energy eigenvalue, and the series extrapolation would no 
longer be valid beyond that point. 

IV. T-Jz MODEL 

We also study a simplified version of the t-J model, the t-Jz model, which has an Ising 
rather than a Heisenberg spin interaction 

H = J. E(^f ^1 - ^^^^^.O - ^ E ncU,,. + ci.c,.)P (21) 

Here we naturally take the first term in H as the unperturbed Hamiltonian, and the second 
term as the perturbation, and perform an expansion in t/J^. This expansion only converges 
well up to about t/J^^I. To improve the convergence in the large t/J^ region, we can 
include a staggered z-field of strength r as in the t-J model, and perform the following 
separation, 

H/J, = Ho + xV (22) 

where 

^0 = E(^r^| - \n,n,)+rY^{-iyS^ (23) 

(ij) i 

V = -ji: ^(4c,.. + 4c,.)P - r E(-l)*^r (24) 

We then take Hq as the unperturbed Hamiltonian and V as the perturbation, and perform 
an expansion in x for fixed values of t/Jz and r. Here again we need to extrapolate the 
series to x = 1. 

The series for the one-hole state have been computed up to order 21 (that is, to (t/Jz)'^^ 
for the Hamiltonian in Eq. (^) or to order for the Hamiltonian in Eq. (p^)) involving a 



^There is a discrepancy between Eder's results and ours, which is due to his omission of the term 
— l/AniUj in the Hamiltonian. 
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list of 23546 clusters with up to 11 sites, while the series for the two-hole bound state have 
been computed up to order 18, involving a list of 23546 clusters up to 11 sites. The final 
series in t/J^ for the spectrum of the one-hole state and two-hole bound state are listed in 
Table |IT|. 

Fig. ^ shows the one-hole spectrum for some particular values of t/ Jz- Here we can see 
that the shape of the dispersion curve is quite different from that for the t-J model: the 
lowest energy occurs at (0, 0) rather than (7r/2, 7r/2) as in the t-J model, and the bandwidth 
here is about 10 times smaller than that for the t-J model. The minimum one-hole energy 
as a function of Jz/t is shown in Fig. |TD|. In the interval 0.2 < Jz/t < 1, the minimum 



one-hole energy can be very well fitted |^1],§ as 

E^Vt = -3.60 + 2.9(J,/t)^/^ + Jz/t (25) 
where the linear term comes from the —l/Ariinj term in the Hamiltonian. This is in good 



agreement with previous estimates |41 



The two-hole dispersion is shown in Fig. |Tl|. We note that the energy at (0, 0) is exactly 
7Jz/2, independent of t, and the energy along path (vr, 0) (p-wave pair) to (vr, tt) (ci-wave 
pair) is very flat. 

The S-, p- and d-wave binding energies defined by 



Eb(k) = E^^{k) - 2E^^{k = 0) (26) 



are shown in Figure |12[ We can see from this figure that s-wave binding disappears around 
t/Jz — 0.32. The d-wave binding energies are slightly greater than for p-wave when t/Jz ^ 
0.83, but when t/Jz ^ 0.83 , the p-wave pair has a greater binding energy than the d-wave, 
similarly to that found for the full t-J model. Thus the t-Jz model also shows a crossover 
between rf-wave and p-wave bound states, occurring even earlier than in the full t-J model. 
Again, the earlier variational calculations of Eder showed very similar results (Fig. 12 of 
Ref. H). 

The solid points in Figure ^ are the results of a reduced basis diagonalization study 



on clusters of up to 50 sites by Riera and Dagotto [^. The trend of these results agrees 
extremely well with the series estimates, although their estimates lie systematically a little 
lower. They do not discuss the symmetry of the bound state They find the binding of 
the 2-hole state disappears at Jz/t ~ 0.183, or t/Jz — 5.5, which agrees very well with our 
results. 



V. CONCLUSIONS 

The t-J model on the square lattice has been studied by many workers over the last 
decade, using a variety of techniques. These methods are all approximate and have uncertain 
accuracy. Independent confirmation from a variety of approaches increases one's confidence 
in the results. 

We have used the linked-cluster series expansion method to study one and two hole 
properties of the half-filled t-J and t-Jz models in 2-dimensions at T = 0. To our knowledge 
this is the first time this has been done. Our results for 1-hole and 2-hole states can be 
summarized as follows: 
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The 1-hole dispersion curve has been computed throughout the Brillouin zone, and 
agrees in shape with previous calculations. 



The minimum for 1-hole dispersion occurs at 
in agreement with previous work. 



TT 



/2,7r/2) and equivalent points, again 



The bandwidth W for the 1-hole dispersion curve of the t-J model has been calculated. 



At small J/ 1 it agrees with the spin- wave calculations of Martinez and Horsch |^T|] and 
Liu and Manousakis [l^, but beyond J/t ~ 1 our results lie substantially higher than 
theirs. We get a peak value of W/t ^ 1.15 at J/t ~ 1.25, rather than W/t ~ 1 at 
J/t c:^ 0.8 as given in previous work. 

We have computed the dispersion curve for 2-hole bound states. To our knowledge 
this has not been well calculated previously. 

In the limit t/J = of the t-J model, we can compare our results for the one-hole 
and two-hole energy {E^'' and E^^), the binding energy (-Eb), and the bandwidth 



for the one- hole state {W) with previous spin- wave [^], Monte Carlo and other 
calculations. 





This Work 


Previous works 


E^^/J 
E^^/J 

W 


2.164(2) 
3.798(2) 
-0.529(4) 
3.26(t/J)2 


2.155(SW), 2.193(MC) 
3.79436(SW), 3.801(8)(MC) 
-0.51(SW), -0.58(2)(MC) 

2{t/jf 0, 2.6(t/j)2 g 



We can see that our results agree very well with other results, except for W , where we 
find a significantly larger value for the leading coefficient. 

At large t / J, we find a crossover in binding energy between the d-wave and p-wave 
2-hole states, for both the t-J and t-J^ models. This phenomenon is likely to arouse 
some controversy: it certainly needs to be checked by other methods. As supporting 
evidence, we note the following points: 

i) The trend of the series results for the t-Jz model agrees very well with numerical 



calculations of Riera and Dagotto ^1| ; 

ii) The phenomenon would be in accord with the Kohn-Luttinger effect as predicted 
by Baranov, Kagan and Chubukov p6| , p7| , and confirmed at low densities p9| , |38| , ^ . 
Chubukov pointed out that the same effect should occur in an arbitrarily dense Fermi 
liquid. 

iii) Very much the same phenomenon was observed previously in the variational cal- 
culations of Eder [|T4|. The concordance between his results and ours would seem to 
indicate that his variational ansatz is a good one. 

iv) A similar enhancement of pairing in the odd-parity singlet channel, analogous to 
our p-wave state, was also noticed for the Hubbard model by Scalettar, Singh and 
Zhang . 

Another possibility, however, is that another state in the ci-wave sector may cross over 
the one we have tracked, and become the lowest energy state at large t/J. Such a 
phenomenon would not be detected by our series analysis. 
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The formation of a two-hole bound state is not of itself sufficient to prove the existence 
of either a superconducting condensate, or of "phase separation" between filled and empty 
sites ; but one would certainly presume it to be a necessary precondition for either of these 
phenomena to occur. Our results indicate that the (i-wave binding energy vanishes at around 
t/J ~ 2.5. The p-wave binding energy may continue a little longer, but it is certainly very 
weak, and our error bars are not small enough to allow any definite conclusion on this point. 
Boninsegni and Manousakis and Poilblanc et al. |0] found the binding to disappear at 
around t/J ~ 3.7. 

The question of phase separation is still under debate: see the review by Dagotto [|I|. 
At small t/J phase separation is known to occur. Emery, Kivelson and Lin |^ suggested 
that in fact phase separation occurs at all couplings in the neighborhood of half-filling. 
Subsequently, the majority of studies p|j39|] concluded that at intermediate coupling (0.3 < 
J/t < 1) there is a region where binding of pairs occurs, but no phase separation. Poiblanc 
and Kohno |^ put the phase separation boundary at about t/J ~ 1.3. But Hellberg 
and Manousakis |39| have recently challenged this scenario, arguing that the criterion of 
a vanishing inverse compressibility used in previous studies was unreliable. They use a 
Maxwell construction technique, from which they conclude that indeed phase separation 
occurs at all couplings in accordance with the original hypothesis of Emery et al. [^]. Our 
results unfortunately have no direct bearing on this debate. To explore these questions one 
must treat the model away from half-filling, which is not possible using our present series 
technique, because it cannot handle an unperturbed ground state which is "doped" away 
from half-filling, and is therefore inhomogeneous. 

The presence or absence of a superconducting condensate is also a matter of debate 
imi. One recent study |^ found no evidence of (i-wave superconducting correlations in the 
physical parameter range t/J > 2. The simple t-J model is beginning to seem a little 
doubtful ClS db theory of the high-Tc superconductors. 

There are a number of other significant physical quantities, such as spectral functions, 
spectral weights, and magnetic susceptibilities which are the subject of ongoing work. 
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FIGURES 

FIG. 1. (a) The open (full) points represent the up (down) sublattices, respectively. The 
primitive translation vectors, {ex,ey) and {ex',ey'), are shown by the arrows, (b) The reciprocal 
space of these lattices. The Brillouin zones for {ex,ey) and {ex',ey') are — vr/a < kx,ky < n/a and 
— 7r/(\/2a) < kx', kyi < 7r/(\/2a), respectively, where a is the lattice spacing. 

FIG. 2. The extrapolation of the j;-series for one-hole energy E^^/ Jz for y = 0.5 at momentum 
k = (vr/2, 7r/2) and several different r values. Several different integrated differential approximants 
to the series are shown. The final estimate at x = 1 is shown as a cross with error bar. 

FIG. 3. Plot of one-hole dispersion A/ J (where A = E^^{k)- E^^{0, 0)) along the path TMXT 
in the Brillouin zone (see Fig. ||(b)) for the system with coupling ratios t/J =1/10, 1/5, 1/3, 1/2, 
2/3 and 1. 

FIG. 4. Similar to Fig. ^ but with the hole dispersion curve plotted along the path rS'M'F in 
the Brillouin zone. 

FIG. 5. The normalized bandwidth W/t = [E^'^iO, 0) - E^^{Tr/2,7r/2)]/t as a function of J/t. 
The solid square points with error bars connected by a solid line are our estimates from x-series, 
while the dashed line are our results in Eq.(p^ obtained from the two- variable expansions. The 



crosses are the estimates of Martinez and Horsch |31], the solid circles are the estimates of Liu and 
Manousakis 



FIG. 6. (a) The open (full) points represent the up (down) sublattices, respectively, while 
the crosses represent the centre points of bonds connecting two adjacent sites. The primitive 
translation vectors, (e^,ey) and {ex',eyi), are shown by the arrows, (b) The reciprocal space of 
these lattices. The Brillouin zones for {ex,ey) and {ex',eyi) are — \/2vr/a < kx,ky < \/27r/a and 
— 7r/(-v/2a) < kx',kyi < 7r/(\/2a), respectively. 

FIG. 7. Plot of 2-hole dispersion £^^'^(k)/J along high-symmetry cuts through the Brillouin 
zone for the system with coupling ratios J/t = 10, 5, 3, 2, 3/2. 



FIG. 8. Plot of the estimated binding energy E\^/ J for s-, p- and d-wave pair states (marked 
by s, p and d, respectively) as functions of t/J. The crosses with error bars connected by dashed 
lines are our estimates from x-series, while the solid lines in the small t/J region are our results 
in Eq.(pO|) obtained from the two- variable expansions. The three solid circles with error bars are 
8x8 Green's function Monte Carlo results taken from Ref. ||l^, the five solid triangles are the 
results of 32-site finite lattice calculations [46,pH|, while the two open circles are the extrapolation 



to bulk limit from finite size calculations [pq 
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FIG. 9. Plot of the one-hole dispersion A/J^ (where A = ^^''(k) - ^^''(0,0)) along the path 
TMXT in the Brillouin zone (see Fig. |^(b)) for the t-Jz model with coupling ratios t/ Jz=5/3, 5/4, 
1 and 0.8. 

FIG. 10. Plot of the one-hole energy at k = for the t — Jz model as a function of 
{Jz/tf/^/[l + {Jz/tf/^]. The solid curves in the large {Jz/t)^^^/[l + (Jz/tf^^] region are the 
integrated differential approximants to the t/Jz series, the points with error bars are estimates 
from x-series, while the solid curve in the small (Jz/t)"^^^ /[i- + (Jz/t)"^^^] region is the least square 
fit E^^/t = -3.60 + 2.9(J2/t)2/3 + J^/t to the results in the interval 0.2 < Jz/t < 1. 

FIG. 11. Plot of the 2-hole dispersion E^^{]i.)/ Jz along high-symmetry cuts through the 
Brillouin zone for the t-Jz model with coupling ratios t/Jz =0.2, 0.4, 0.6, 0.8 and 1. 

FIG. 12. Plot of the binding energy E\y/J as a function of t/Jz for the t-Jz model. The 
solid lines are the different integrated differential approximants to the expansion in t/Jz, while the 
points with error bars connected by dashed lines are the estimates from the expansion in x. The 
solid points are the results of a finite size calculations . 
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TABLES 



TABLE I. Series coefficients for tlie one-fiole energy E^^/J^ = J2k^kX^ ^'^^ V — ^-^ 
k = (7r/2, 7r/2) for the t-J model. Nonzero coefficients Cj up to order A; = 11 for r = 0, 1, 2, 4, 8 are 
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TABLE II. Scries coefficients for tlie one-liole and two-hole energy E{kx, ky) /Jz =J2knm '^k,n,m 
x^[cos{nkxA) cos{mkyA) + cos{mkxA) cos{nkyA)]/2 of the t-J model, where A = a for the one-hole 
case and A = aj\f2 for the two-hole case. Nonzero cocfhcicnts (ik.n.m up to order k = \\ for y = 0.5 
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1.5.>1781787x 10 


-2 


( 3, 


1, 


1) 


-3.33.33.33.33.3x10- 


1 


(4, 


1, 


0) 


-2.67.3809524x10" 


1 


(11, 


5, 


4) 


-7.959003190x10" 


~y 


(8, 


4, 


1) 


2.212704012x10" 


-3 


( 4, 


1, 


1) 


-1.192857143x10" 


1 


(5, 


1, 


0) 


-4.024977954x10" 


1 


( 4, 


2, 


0) 


2.962962963x10" 


2 


(9, 


4, 


1) 


-8.862200427x10" 


-4 


( 5, 


1, 


1) 


2.904905203x10" 


1 


(6, 


1, 


0) 


-1.. 302727805x10" 


1 


( 5, 


2, 


0) 


-2.629629630x10" 


2 


(10, 


4, 


1) 


-9.978608234x10" 


-2 


( 6, 


1, 


1) 


5.329329186x10" 


1 


(7, 


1, 


0) 


3.4995.36691x10" 


1 


( 6, 


2, 


0) 


-2.375711537x10" 


1 


(11, 


4, 


1) 


1.15.3622881x10" 


-1 


( 7, 


1, 


1) 


-6.271803816x10" 


1 


(8, 


1, 


0) 


4.079418136x10" 


1 


( 7, 


2, 


0) 


2..330090787X10" 


1 


(10, 


5, 


2) 


4.230241698x10" 


-4 


( 8> 


1, 


1) 


-7.522337850x10" 


1 


(9, 


1, 


0) 


5..301735104X10" 


2 


( 8, 


2, 


0) 


3.039557383x10" 


1 


(11, 


5, 


2) 


-1.750974021x10" 


-4 


( 9, 


1, 


1) 


5..347837206X10" 


1 


(10, 


1, 


0) 


-1.026434348 




( 9, 


2, 


0) 


-3.529882617x10" 


1 


(8, 


4, 


0) 


1.106352006x10" 


-3 


(10, 


1, 


1) 


1.146091910 




(11, 


1, 


0) 


-7.412817539x10" 


1 


(10, 


2, 


0) 


-6.813379022x10" 


1 


(9, 


4, 


0) 


-7.128118497x10" 


-4 


(11, 


1, 


1) 


-4.863394606x10" 


1 


(4, 


2, 


1) 


5.925925926x10" 


2 


(11, 


2, 


0) 


9.790814131x10" 


1 


(10, 


4, 


0) 


-5.5388.38348x10" 


-2 


( 4, 


2, 


2) 


1.481481481x10" 


2 


(5, 


2, 


1) 


-2.629629630x10" 


2 


( 6, 


3, 


1) 


8.747795414x10" 


3 


(11, 


4, 


0) 


8.043954590x10" 


-2 


( 6, 


2, 


2) 


-5.533170820x10" 


2 


(6, 


2, 


1) 


-8.610412317x10" 


-2 


( 7, 


3, 


1; 


-5.867614638x10" 


3 


(10, 


5, 


1) 


5.287802122x10" 


-4 


( 7, 


2, 


2) 


1.147526740x10" 


2 


(7, 


2, 


1) 


5.769980192x10" 


-2 


( 8. 


3, 


1; 


-1.459267389x10" 


1 


(11, 


5, 


1) 


-2.626461032x10" 


-4 


( 8, 


2, 


2) 


3.484612337x10" 


2 


(8, 


2, 


1) 


2.118965896x10" 


-1 


( 9. 


3, 


1; 


1.611775400x10" 


1 


(10, 


5, 


0) 


3.172681273x10" 


-4 


( 9> 


2, 


2) 


-8.489781072x10" 


2 


(9, 


2, 


1) 


-9.184491129x10" 


-2 


(10, 


3, 


1; 


1.623787831x10" 


1 


(11, 


5, 


0) 


-9.550737020x10" 


-5 


(10, 


2, 


2) 


-8.365174384x10" 


2 


(10, 


2, 


1) 


-3.638022235x10" 


1 


(11, 


3, 


1) 


-3.687197934x10" 


1 












(11, 


2, 


2) 


3.945274143x10" 


1 


(11, 


2, 


1) 


-5.611998341x10" 


3 


( 8, 


4, 


2) 


1.475136008x10" 


3 
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TABLE III. Series coefficients for the one-hole and two-hole energy E{kx, ky)/ Jz = J2k,n,m 
cik,n,m {t/ Jz)^[cos{nkxA) cos{'mkyA) + cos{'mkxA) cos{nkyA)]/2 for the Hamiltonian in Eq. (^Tj) of 
the t-Jz model, where A = a for the one- hole case and A = a/V2 for the two-hole case. Nonzero 



coefficients ak^n,m up 


to order k = 


20 for one-hole 


case OT k = 


18 for two-hole case are listed. 


(k,n,m) afc,,i,m 


1 (k,n,m) 


^k.n,m 1 


(k,n,m) 


ak,n,m 1 (k,n,m) afe,n,m 



one hole energy E'-^(kx , ky)/ Jz 



(0, 


0, 


0) 


2.000000000 




(20, 


0, 


0) 


3.269830958x10* 


(8, 


1, 


1) 


1.107677837 


(20, 3, 


1) 


1.363566623 




(2, 


0, 


0) 


-2.666666667 




(10, 


2, 


0) 


-9.187011686x10"^ 


(10, 


1, 


1) 


-5.258263841 


(12, 2, 


2) 


2.213826787x10" 


2 


(4, 


0, 


0) 


2.607407407 




(12, 


2, 


0) 


1.237459977x10"! 


(12, 


1, 


1) 


2.361854957x10! 


(14, 2, 


2) 


-2.801396443x10" 


1 


(6, 


0, 


0) 


-5.667818930 




(14, 


2, 


0) 


-1.042288265 


(14, 


1, 


1) 


-1.041949744x10^ 


(16, 2, 


2) 


2.293697931 




(8, 


0, 


0) 


1.579276071x10^ 




(16, 


2, 


0) 


7.139692772 


(16, 


1, 


1) 


4.570753372 xlO^ 


(18, 2, 


2) 


-1.550982309x10! 




(10, 


0, 


0) 


-4.977461910x10^ 




(18, 


2, 


0) 


-4.367374242 xlO^ 


(18, 


1, 


1) 


-2.003192131 xlO^ 


(20, 2, 


2) 


9.435217647x10! 




(12, 


0, 


0) 


1. 689229038 xlO^ 




(20, 


2, 


0) 


2.492266705 xlO^ 


(20, 


1, 


1) 


8.787977137x10^ 


(20, 4, 


2) 


4.053738087x10" 


5 


(14, 


0, 


0) 


-6.023423934x10^ 




(18, 


4, 


0) 


-1.275445883x10"^ 


(14, 


3, 


1) 


-1.112602736x10"* 


(18, 3, 


3) 


-4.160122546x10" 


3 


(16, 


0, 


0) 


2.225264150x10^ 




(20, 


4, 


0) 


1.432818794x10"* 


(16, 


3, 


1) 


6.891013162x10"^ 


(20, 3, 


3) 


8.089745849x10" 


2 


(18, 


0, 


0) 


-8.447437043 xlO^ 




(6, 


1, 


1) 


-1.896296296x10"! 


(18, 


3, 


1) 


-1.208165780x10"! 




























two-hole energy 




kx ; k 


y)l 












(0, 


0, 


0) 


3.500000000 




(2, 


1, 


0) 


2.666666667 


(10, 


3, 


1) 


-3.833568660 


(10, 


4, 


0) 


-4.187749256x10" 


1 


(2, 


0, 


0) 


-4.000000000 




(4, 


1, 


0) 


-2.607407407 


(12, 


3, 


1) 


3.749428458x10! 


(12, 


4, 


0) 


-4.625682644 




(4, 


0, 


0) 


2.311111111 




(6, 


1, 


0) 


-9.124044680x10"! 


(14, 


3, 


1) 


-1.763568623x102 


(14, 


4, 


0) 


3.879009243x10! 




(6, 


0, 


0) 


-3.578130511x10" 


1 


(8, 


1, 


0) 


2.709462265 xio! 


(16, 


3, 


1) 


5.362296757x10^ 


(16, 


4, 


0) 


-1.755530196x10^ 




(8, 


0, 


0) 


-1.447791909x101 




(10, 


1, 


0) 


-1.454354060x10^ 


(18, 


3, 


1) 


-4.338040480x10^ 


(18, 


4, 


0) 


5.141129944x10^ 




(10, 


0, 


0) 


8.449051521 xlO^ 




(12, 


1, 


0) 


4.973154521 xlO^ 


(8, 


4, 


2) 


3.776348180x10"! 


(10, 


5, 


1) 


5.414709373x10" 


1 


(12, 


0, 


0) 


-3.205718341 xlO^ 




(14, 


1, 


0) 


-8.044807264x10^ 


(10, 


4, 


2) 


-7.694903208x10"! 


(12, 


5, 


1) 


-6.844512414x10" 


1 


(14, 


0, 


0) 


7.683946158 xlO^ 




(16, 


1, 


0) 


-3.592648436 xlO^ 


(12, 


4, 


2) 


-3.834200865 


(14, 


5, 


1) 


-1.368955321x10! 




(16, 


0, 


0) 


1.864995221 xlO^ 




(18, 


1, 


0) 


4.011882196x10* 


(14, 


4, 


2) 


4.301073460x10! 


(16, 


5, 


1) 


1.169642314x10^ 




(18, 


0, 


0) 


-1.447204343x10* 




(4, 


2, 


1) 


9.481481481x10"! 


(16, 


4, 


2) 


-2.363456039 xlO^ 


(18, 


5, 


1) 


-5.467419040 xlO^ 




(2, 


1, 


1) 


1.333333333 




(6, 


2, 


1) 


-2.072004703 


(18, 


4, 


2) 


9.094821992 xlO^ 


(12, 


6, 


2) 


5.024245295x10" 


1 


(4, 


1, 


1) 


-1.362962963 




(8, 


2, 


1) 


1.071796284 


(10, 


5, 


3) 


2.707354687x10"! 


(14, 


6, 


2) 


-5.607740036x10" 


1 


(6, 


1, 


1) 


3.117601411 




(10, 


2, 


1) 


1.934131152x10! 


(12, 


5, 


3) 


-6.099856791x10"! 


(16, 


6, 


2) 


-1.688448287x10! 




(8, 


1, 


1) 


-1.157205899x10^ 




(12, 


2, 


1) 


-1.610227509x10^ 


(14, 


5, 


3) 


-3.591949866 


(18, 


6, 


2) 


1.530376091 xlO^ 




(10, 


1, 


1) 


3.724437583 xlO^ 




(14, 


2, 


1) 


9.100666397x10^ 


(16, 


5, 


3) 


4.680082984x10! 


(14, 


7, 


3) 


4.576926174x10" 


1 


(12, 


1, 


1) 


-4.436033620x10^ 




(16, 


2, 


1) 


-4.053430224x10^ 


(18, 


5, 


3) 


-2.944841481 xlO^ 


(16, 


7, 


3) 


-4.694339733x10" 


1 


(14, 


1, 


1) 


-5.334708723 xlO^ 




(18, 


2, 


1) 


1.279544679x10* 


(12, 


6, 


4) 


2.009698118x10"! 


(18, 


7, 


3) 


-1.886986003x10! 




(16, 


1, 


1) 


5.421019234x10^ 




(6, 


3, 


2) 


3.732392710x10"! 


(14, 


6, 


4) 


-4.925030503x10"! 


(16, 


8, 


4) 


4.117367040x10" 


1 


(18, 


1, 


1) 


-3.093983511x10* 




(8, 


3, 


2) 


-7.739781828x10"! 


(16, 


6, 


4) 


-3.212472559 


(18, 


8, 


4) 


-4.052076547x10" 


1 


(4, 


2, 


2) 


2.370370370x10" 


1 


(10, 


3, 


2) 


-1.053007756 


(18, 


6, 


4) 


4.903204193x10! 


(18, 


9, 


5) 


3.670192948x10" 


1 


(6, 


2, 


2) 


-5.237389771x10" 


1 


(12, 


3, 


2) 


1.967483048x10! 


(14, 


7, 


5) 


1.525642058x10"! 


(10, 


5, 


0) 


3.248825624x10" 


1 


(8, 


2, 


2) 


1.010441376 




(14, 


3, 


2) 


-1.305817944x10^ 


(16, 


7, 


5) 


-4.034853659x10"! 


(12, 


5, 


0) 


-6.899424873x10" 


1 


(10, 


2, 


2) 


-1.494008026 




(16, 


3, 


2) 


6.744405152 xlO^ 


(18, 


7, 


5) 


-2.770916581 


(14, 


5, 


0) 


-6.463536863 




(12, 


2, 


2) 


1.846909265x10" 


1 


(18, 


3, 


2) 


-2.972249754x10^ 


(16, 


8, 


6) 


1.176390583x10"! 


(16, 


5, 


0) 


6.567477615x10! 




(14, 


2, 


2) 


1. 042676508 xlOi 




(8, 


4, 


3) 


1.888174090x10"! 


(18, 


8, 


6) 


-3.342928209x10"! 


(18, 


5, 


0) 


-3.388841937x10^ 




(16, 


2, 


2) 


-1.111625374x10^ 




(10, 


4, 


3) 


-2.886375899x10"! 


(18, 


9, 


7) 


9.175482371 xl0"2 


(12, 


6, 


1) 


6.698993726x10" 


1 


(18, 


2, 


2) 


1.249271047x10^ 




(12, 


4, 


3) 


-1.577708121 


(6, 


3, 


0) 


3.732392710x10"! 


(14, 


6, 


1) 


-1.079421253 




(6, 


3, 


3) 


9.330981775x10" 


2 


(14, 


4, 


3) 


1.855760450x10! 


(8, 


3, 


0) 


-9.806190570x10"! 


(16, 


6, 


1) 


-2.116718725x10! 




(8, 


3, 


3) 


-2.869718840x10" 


1 


(16, 


4, 


3) 


-1.206715261 xlO^ 


(10, 


3, 


0) 


-7.606578761x10"! 


(18, 


6, 


1) 


2.016971553x10^ 




(10, 


3, 


3) 


7.302649849x10" 


1 


(18, 


4, 


3) 


6.228213110x10^ 


(12, 


3, 


0) 


1.900360917x10! 


(14, 


7, 


2) 


6.538465963x10" 


1 


(12, 


3, 


3) 


6.474454323x10" 


2 


(10, 


5, 


4) 


1.082941875x10"! 


(14, 


3, 


0) 


-1.096944607x10^ 


(16, 


7, 


2) 


-6.381513372x10" 


1 


(14, 


3, 


3) 


-1.507976852x10^ 




(12, 


5, 


4) 


-7.182655253xl0"2 


(16, 


3, 


0) 


4.525126033 xlO^ 


(18, 


7, 


2) 


-2.878600133x10! 




(16, 


3, 


3) 


1.112739551 xlO^ 




(14, 


5, 


4) 


-1.662709009 


(18, 


3, 


0) 


-1.604127694x10^ 


(16, 


8, 


3) 


6.176050560x10" 


1 


(18, 


3, 


3) 


-5.501017323x10^ 




(16, 


5, 


4) 


1.713732791x10! 


(8, 


4, 


1) 


5.664522270x10"! 


(18, 


8, 


3) 


-1.660697208x10" 


1 


(8, 


4, 


4) 


4.720435225x10" 


2 


(18, 


5, 


4) 


-1.134996640x10^ 


(10, 


4, 


1) 


-1.358637239 


(18, 


9, 


4) 


5.709189031x10" 


1 


(10, 


4, 


4) 


-1.753576976x10" 


1 


(12, 


6, 


5) 


6.698993726x10"^ 


(12, 


4, 


1) 


-5.476288758 


(12, 


6, 


0) 


3.349496863x10" 


1 


(12, 


4, 


4) 


5.665451332x10" 


1 


(14, 


6, 


5) 


2.348849413x10"^ 


(14, 


4, 


1) 


6.555818772x10! 


(14, 


6, 


0) 


-2.248529752x10" 


1 


(14, 


4, 


4) 


7.031489323x10" 


1 


(16, 


6, 


5) 


-1.498021733 


(16, 


4, 


1) 


-3.561148091 xlO^ 


(16, 


6, 


0) 


-1.304144759x10! 




(16, 


4, 


4) 


-2.234394856 xlO^ 




(18, 


6, 


5) 


1.514875202x10! 


(18, 


4, 


1) 


1.327749407x10^ 


(18, 


6, 


0) 


1.069591051 xlO^ 




(18, 


4, 


4) 


1.696295810x10^ 




(14, 


7, 


6) 


4.358977309x10"^ 


(10, 


5, 


2) 


4.331767498x10"! 


(14, 


7, 


1) 


7.628210290x10" 


1 


(10, 


5, 


5) 


2.707354687x10" 


2 


(16, 


7, 


6) 


6.397044909 xl0"2 


(12, 


5, 


2) 


-8.287589771x10"! 


(16, 


7, 


1) 


-1.648715187x10" 


1 


(12, 


5, 


5) 


-1.145505796x10" 


1 


(18, 


7, 


6) 


-1.267992012 


(14, 


5, 


2) 


-8.117281873 


(18, 


7, 


1) 


-3.918669310x10! 




(14, 


5, 


5) 


4.544498080x10" 


1 


(16, 


8, 


7) 


2.940976457x10"^ 


(16, 


5, 


2) 


8.399862105x10! 


(16, 


8, 


2) 


8.234734079x10" 


1 


(16, 


5, 


5) 


1.078497914 




(18, 


8, 


7) 


7.770975088 xl0"2 


(18, 


5, 


2) 


-4.609565870x10^ 


(18, 


8, 


2) 


1.073653528x10" 


1 


(18, 


5, 


5) 


-2.724985485 xlO^ 




(18, 


9, 


8) 


2.038996082 xl0"2 


(12, 


6, 


3) 


3.349496863x10"! 


(18, 


9, 


3) 


8.563783546x10" 


1 


(12, 


6, 


6) 


1.674748432x10" 


2 


(4, 


2, 


0) 


4.740740741x10"! 


(14, 


6, 


3) 


-4.223690648x10"! 


(14, 


7, 


0) 


4.358977309x10" 


1 


(14, 


6, 


6) 


-7.829101098x10" 


2 


(6, 


2, 


0) 


-6.512874780x10"! 


(16, 


6, 


3) 


-9.671323639 


(16, 


7, 


0) 


-4.703472705x10" 


1 


(16, 


6, 


6) 


3.741008913x10" 


1 


(8, 


2, 


0) 


-1.556811959 


(18, 


6, 


3) 


9.548744493x10! 


(18, 


7, 


0) 


-1.979896930x10! 




(18, 


6, 


6) 


1.267008305 




(10, 


2, 


0) 


1.207544510x10! 


(14, 


7, 


4) 


2.615386385x10"! 


(16, 


8, 


1) 


1.029341760 
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(14, 7, 7) 1.089744327x10-2 
(16, 7, 7) -5.528545914x10-2 
(18, 7, 7) 3.138641071x10-1 

(16.8.8) 7.352441142x10-'^ 
(18, 8, 8) -4.001373177x10-2 

(18.9.9) 5.097490206x10-3 



(12, 2, 0) -3.196277031 xlO^ 
(14, 2, 0) -5.271684834x10^ 
(16, 2, 0) 1.018254901x103 
(18, 2, 0) -5.459613629x103 
( 6,3,1) 5.598589065x10-1 
( 8, 3, 1) -9.918360628x10-1 



(16, 7, 4) -1.332433821x10-1 
(18, 7, 4) -1.033273379x101 
(16,8,5) 2.058683520x10-1 

(18.8.5) 6.203359335x10-2 

(18.9.6) 1.631196866x10-1 
( 8,4,0) 2.832261135x10-1 



(18, 8, 1) -4.479185205x10-1 
(18, 9, 2) 1.141837806 
(16,8,0) 5.146708800x10-1 
(18,8,0) 2.583076501x10-1 
(18, 9, 1) 1.284567532 
(18,9,0) 7.136486288x10-1 
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